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Abstract 

We explore viable scenarios for parametric resonant amplification 
of electroweak (EW) gauge fields and Chern-Simons number during 
preheating, leading to baryogenesis at the electroweak (EW) scale. 
In this class of scenarios time-dependent classical EW gauge fields, 
essentially spatially-homogeneous on the horizon scales, carry Chern- 
Simons number which can be amplified by parametric resonance up 
to magnitudes at which unsuppressed topological transitions in the 
Higgs sector become possible. Baryon number non-conservation asso- 
ciated with the gauge sector and the highly non-equilibrium nature of 
preheating allow for efficient baryogenesis. The requisite large CP vi- 
olation can arise either from the time-dependence of a slowly varying 
Higgs field (spontaneous baryogenesis), or from a resonant amplifica- 
tion of CP violation induced in the gauge sector through loops. We 
identify several CP violating operators in the Standard Model and its 
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minimal extensions that can facilitate efficient baryogenesis at pre- 
heating, and show how to overcome would-be exponential suppression 
of baryogenesis associated with tunneling barriers. 

1 Introduction 

Parametric resonance coupling of the oscillating inflaton to Standard Model 
fields [m offers new dynamical mechanisms for various early-universe phe- 
nomena. One interesting scenario 0, ^, Q adjusts inflaton parameters^ so 
that reheating does not heat the universe to a temperature above the EW 
cross-over temperature Tc. This means that sphaleron transitions are frozen 
out and that B+L created before reheating will not be washed out. 

The early attempts ||, to create B+L solely through Standard Model 
effects suffered from the fact that if the only source of CP violation is in 
the CKM matrix, it would lead to far too small a value of B+L, and from 
the need for a first-order EW phase transition to drive out-of-equilibrium 
effects. This first-order transition seems to be ruled out by current lower 
limits on the Higgs mass of around 110 GeV. So it might appear that EW 
baryogenesis is unattractive. However, with the suggestion of EW preheating 
and the accompanying parametric resonances, it has become interesting to 
look at EW baryogenesis in a different way. 

This was done in a recent paper 0, where two of us showed that EW para- 
metric resonance, with the Higgs field oscillating because of its coupling to the 
oscillating inflaton, could amplify CP- violating seed values of the EW gauge 
potential in an ansatz where this gauge potential was spatially-homogeneous 
(out to the horizon) but time- dependent. The CP violation was manifested 
as a homogeneous classical "condensate" of Chern-Simons number or equiv- 
alently of B+L (through the EW anomaly). The scenario of P] had several 
shortcomings: 

1. CP violation only occurred because of CP- violating initial values for 
the gauge potential, not because of any explicit CP violation in the 
gauge equations of motion. 



^The consistency of low-scale inflation with primordial density perturbations has been 
shown by German 
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2. The Higgs field oscillations were considered as given, so that no (classi- 
cal) back-reaction of the gauge fields on the Higgs field was considered. 

3. No plausible scenario was suggested for the permanent conversion of 
Chern-Simons number to actual baryons and leptons; once parametric 
resonance driving ended, it would be possible for the homogeneous 
Chern-Simons condensate to dissipate. Actual formation of baryons 
and leptons requires the development of spatial inhomogeneities, rather 
similar to sphalerons; one might term this process sphalerization. 

In the present work we extend the considerations of p[| to deal, at least 
partially, with these three shortcomings. 

First, we introduce explicit CP violation]^ in the EW gauge plus Higgs 
equations of motion. This CP violation comes from coupling of CP-violating 
effects to the EW gauge fields through loops, typically quark loops. We 
explore two cases; in one, there is out-of-equilibrium strong CP violation as 
evidenced by an rj' field slowly (on EW time scales) rolling in a standard 
QCD potential with zero 6 angle (so that there is no strong CP violation 
when the rj' field is in equilibrium). In the other, we invoke spontaneously- 
broken CP in multi-Higgs models 0], again, out of equilibrium.^ Note that 
it is not easy to rule out the CP-violating operators we use by present-day 
experimental results on CP violation, since in our scenario the CP-violating 
effects are far from their equilibrium values. 

Second, we study the combined gauge field-Higgs classical equations, 
extending the previous ansatz to include a spatially-homogeneous time- 
dependent Higgs field. The classical backreaction of the gauge potential 
on the Higgs field can drive the Higgs field into oscillating through zero VEV 
rather than staying near a broken-symmetry minimum. This is important 
for sphalerization. 

Third, in the final section of the paper we make some remarks about the 
dynamics of sphalerization. This is a truly difficult dynamical problem, re- 

^The CP-violating terms used here are similar to those invoked in spontaneous baryoge- 
nesis but represent couplings to the gauge sector; see [D and references therein. We do not 
require a first-order phase transition, as commonly invoked in spontaneous baryogenesis. 
For a version of spontaneous baryogenesis in the present context, see Sec. |^below. 

■^These multi-Higgs models are not realistic for CP violation in the K — K system, but 
they could still play a central role in baryogenesis; note that in general our scenarios for 
baryogenesis involve CP violation which is far from the equilibrium values seen today. 
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quiring extensive numerical investigation, and we barely scratch the surface 
here. Aside from brute-force computation, not attempted here, there are sev- 
eral avenues to explore for guidance, including simple effective-temperature 
arguments previously used in connection with EW preheating [|l], ^ ^; nu- 
merical studies of the 1-1-1- dimensional Abelian Higgs model ^ |l^; and 
approximate but useful tools for an analysis based on tools developed for 
understanding the possibility of B-l-L violation in high-energy two-particle 
collisions (see, e.g., [|12])- A major obstacle to the analysis is that mere 



production of Chern-Simons number is not enough to produce baryons; the 
Chern-Simons number must be converted into baryons through the medium 
of formation of Higgs winding number A^wind- Using simple topological ar- 
guments, we construct a model for the energy profile of the system which 
illustrates the issues involved. With the simplification of using an effective 
temperature we estimate baryoproduction for an arbitrary time dependence 
of the topological transition rate and show how wash-out — fermionic back- 
reaction leading to dissipation of the newly-created baryons and leptons — 
can be taken into account. Finally we attempt an alternative to effective- 
temperature considerations by extending some approximate techniques used 
long ago to examine the possibility of B-l-L violation in two-particle collisions. 

The upshot of these considerations is that there are at least two possibil- 
ities for conversion of Chern-Simons number to baryons which do not suffer 
from the exponential suppression of tunneling. In the first, various effects 
(such as gauge backreaction on the Higgs field) may cause the Higgs VEV 
v{t) to oscillate through zero during or immediately after preheating, allow- 
ing unsuppressed transitions (the sphaleron mass Ms ^ v/g vanishes). In 
the second, formation of spatial structures on various scales may allow for 
baryon formation to proceed at energies above barrier heights even if v is near 
its usual broken-symmetry value. In this case, there is a relation between the 
size of the structure, the energy of the structure, and the Chern-Simons den- 
sity, which we explore with a relatively crude but qualitatively satisfactory 
model of the topological charge barrier factor as a function of energy and 
size. We argue that as various spatial scales arise and grow during preheat- 
ing and the Chern-Simons condensate becomes inhomogeneous, it is possible 
to have conversion of Chern-Simons number to actual baryons (and leptons, 
which we ignore in this paper) which is not only unsuppressed with regard 
to topological charge tunneling, but is also not exponentially suppressed by 
poor overlap between initial and final states. This is an important point. 
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About a decade ago, it was suggested that EW collisions of two particles 
yielding an A^-particle final state at energies E > Ms large enough to over- 
come the barrier height and with ~ 4:7i/g'^, could lead to unsuppressed 
B+L production (see, e.g., [|12| for a contemporaneous collection of papers 
on the subject). But others (including Banks et al [jl2| and Cornwall |jT2|, [l^) 



pointed out that the poor overlap between the initial and final states led to 
exponential suppression anyhow, with a barrier which was a finite fraction 
of the canonical 't Hooft barrier. The suppression mechanism had actually 
been proposed earlier by Drukier and Nussinov [|l4l to show that production 
of solitons such as the sphaleron in two-particle collisions was exponentially 
suppressed. 

In the present case, the initial state for baryogenesis is not a two-particle 
plane- wave state, but a Chern- Simons condensate which is largely homoge- 
neous and on which various spatial ripples are growing, as suggested in [^]; 
when the gauge potential grows to be of order m/g, the spatial growth rate 
is of order m. We show that this sort of condensate can have good over- 
lap with baryonic final states (as connected via the EW anomaly) at certain 
spatial scales under the same circumstances (generally high enough energy) 
where the topological barrier is gone. These spatial scales, to no one's sur- 
prise, are around at energies near the sphaleron mass. Of course, if the 
Higgs VEV V is rather small, the corresponding spatial scale grows larger 
and the energy scale grows smaller. We show that if this second mechanism 
to avoid suppression is to work when m is near its vacuum value, the origi- 
nal spatially-homogeneous Chern-Simons condensate must have a density of 
order 10~^ — 10~^ units of topological charge in a volume m~^. Such values 
are indeed reached in numerical simulations (this paper and Ref. |^). 

We cannot, with these crude approaches, begin to quantify the number of 
baryons produced. All we can say is that the number of baryons is essentially 
linear in the strength of the CP-violating operators we discuss, and that 
we are not aware of any experimental limits which would lead to strengths 
too small to produce the observed number of baryons, provided that the 
baryogenesis process is not exponentially suppressed and that further baryon 
washout at reheating is not strong. Given the results of this paper, the 
unexplored rate-limiting step for baryogenesis is the growth of perturbations 
at various spatial scales. 

For some simple versions of the models we study it is possible to do some 
approximate analysis of growth of Chern-Simons number in the spatially- 
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homogeneous phase. Generally we find that adding explicit CP violation 
to the gauge equations of motion leads, as expected, to secular (that is, 
not oscillating with the infiaton or Higgs) CP-violating terms in the gauge 
potential and Chern-Simons density. If the explicit CP violation is correlated 
on super-Hubble scales, as might be expected following infiation and which 
we will assume for this paper, the secular average will also be correlated 
on super-Hubble scales even though initial values of the gauge potentials 
might be random from one Hubble domain to the next. Without explicit 
CP violation in the equations of motion, CP violation can only come from 
CP- violating initial values of the EW gauge fields; if (as mentioned in ^) 
these are random from one Hubble volume to the next and on the average 
CP-symmetric, the baryon number averaged over the whole universe will be 
zero and fiuctuations will be too small by the usual factor of l/\fN, where 
is the number of Hubble volumes. (Of course, CP-violating effects before 
preheating can be correlated on super-Hubble scales, so that there can be a 
secular average with no explicit CP violation in the equations of motion.) We 



offer one model, somewhat resembling a Brownian ratchet |]15|, in which the 
secular Chern-Simons average depends not only on the sign of the explicit 
CP violation (assumed to have super-Hubble correlations) but also on the 
sign of initial conditions, which we could assume as random and uncorrelated. 
Nonetheless, this model leads to a sort of secular Chern-Simons average, since 
resonant growth can occur only for one sign of initial-condition parameters; 
otherwise, there is damping. In this model, however, the secular average itself 
wanders slowly on a long time scale and consequently is somewhat inefficient. 
We also offer other models in which the secular Chern-Simons sign depends 
only on the sign of the explicit CP-violating term in the equation of motion; 
these lead to a Chern-Simons condensate of unique sign across the universe. 

We will give a few examples of numerical simulations, making no attempt 
to cover the possible range of initial conditions and models. The examples are 
chosen to illustrate strong resonance and hence strong amplification of CP 
violation; in effect, this means that EW gauge potentials become of order the 
vacuum Higgs VEV v and the Chern-Simons density grows to order m^/Svr^. 
In strong resonance these final values are essentially independent of the initial 
values; only the time needed to reach the final values changes with initial 
conditions. With no resonance there is little or no amplification, and (in our 
scenarios) far too few baryons would be produced. We are in no position 
to say what "final" values of the Chern-Simons condensate are needed to 
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produce today's B+L values, since our understanding of the conversion of 
Chern-Simons number to B+L in this very non-equihbrium process is still 
quite primitive. 

Our general conclusions are that while there are still many unknowns, 
involving both parameters of CP-violating physics and unsolved dynamics, 
we know of nothing which would rule out generation of the observed number 
of baryons in the universe today in this EW preheating scenario. 

Other points of view [0, 0, |^, |T6[ have been expressed concerning the mech- 
anisms at work during EW preheating, including the invocation of approx- 
imate thermalization of long-wavelength gauge boson modes at an effective 
temperature much larger than the eventual reheat temperature, so that EW 
symmetry is effectively restored, and more or less genuine thermalization of 
shorter wavelengths. Without further numerical simulations it is impossible 
for us to say how such ideas compare to what we suggest here, which is based 
on ultimate conversion of a spatially-homogeneous Chern-Simons condensate 
[§] to a condensate more resembling sphalerons. 

There is one numerical lattice simulation of the full d=3+l Higgs- 
gauge system which does not show any secular growth of Chern-Simons num- 
ber, which instead decays to near zero. For technical reasons these authors 
did not include CP-violating terms in the gauge equations of motion. As we 
show here, such terms are important to establish a long-term Chern-Simons 
condensate, which may undergo sphalerization and conversion to baryons. 



2 Spontaneous baryogenesis at preheating 

Biasing the baryon asymmetry through an effective chemical potential can be 
achieved in a model with two Higgs doublets in what is known as spontaneous 
baryogenesis^ Cohen, Kaplan, and Nelson 0, |] proposed that the effective 
T-reversal asymmetry may come from a time dependence in the solution for 
the Higgs field. Their scenario used the variation of the Higgs field inside 
a wall of a bubble formed in a first-order phase transition. A similar effect 
can occur at preheating uniformly in space, on the horizon scales. We will 
discuss this scenario on an example of a two-doublet model. 

The Higgs potential in a general model with two doublets. Hi and 
has a form 

^ For numerical simulations of baryogenesis in two-Higgs models see jl^ . 
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V{H,,H2) = Xi(HlH,-vlf 

+ U{HlH, - vl) + {HlH^ - vl)f 

+ x4{hIhi){hIh2) - {hIh2){hIhi)] 

+ X^lReiHlH^) - V1V2 cos e]' 

+ X6[lm{HlH2)-viV2smC]^ (1) 

At finite temperature, all and Vi receive thermal corrections and depend 
on the temperature. 

During preheating the Higgs fields move along some classical trajectory 

Hi = A(t)e^^^W (2) 

that satisfies the equations of motion 

0. + 3he, + ^ + pr'^ = 0, 
Pi o9i 

■ 2 dV 

Pi + ?>hpi -Oi pi + — = 0, (3) 
opi 

where h is the Hubble constant. The term [dV/dOi) is non- vanishing as long 
as there is a mixing between the Higgs bosons; this term is periodic in 6i. 

In the course of reheating, fermions are created and a thermal equilibrium 
is achieved at some temperature Tr < lOOGeV, low enough to prevent any 
sphaleron transitions. The Higgs fields change from their zero-temperature 
values at the end of infiation to some temperature-dependent VEV: 

at T = 0, pi = Vi] (4) 
atr = T^, pi = Vi{T^). (5) 

At the same time, the phase 9 also changes: 

^(0) ^ em-02{S)) = i, (6) 

^(TJ ^ ^i(Tj-^2(Tj=e(Tj. (7) 
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The time derivative of 6 serves as a chemical potential for the baryon 
number ||^, |T9| because of an effective coupling 



{doe)f, (8) 

that appears in the Lagrangian after the time-dependent phase is eliminated 
from the Yukawa couplings. Here = iritip'-j^'il) + ... is the fermionic current. 

If the gauge fields grow in resonance as described in Ref . 0] , but the Higgs 
fields are out of resonance, the p/ p term in equation can be neglected. The 
equation allows for a slowly varying solution 6{t) that interpolates between 
= and 6 = ^(T^)- Since the Chern-Simon number violating processes 
in the gauge sector are very rapid on the time scale of thermalization, the 
fermions produced during reheating will have time to equilibrate to the min- 
imum of free energy when a thermal distribution is ultimately achieved. The 
effective chemical potential Pg is proportional to 6, and the equilibrium value 
of baryon asymmetry is 

10-° Tl (1^) , (9) 

where is the time of reheating and is the Hubble time at the electroweak 
scale, yt is the top Yukawa coupling, and v(T) is the Higgs VEV. For < 
70GeV, there is no wash-out of the baryon number by thermal sphalerons. 
The required baryon asymmetry, 10~^°, can be achieved in this scenario if 
the reheat time is 10~^ the Hubble time. The reheat time is usually much 
shorter than the Hubble time, and the requisite ratio 10^^ can be achieved 
in a realistic model. 

The difference with the scenario proposed by Cohen, Kaplan and Nel- 
son is that in our case CP violation occurs homogeneously in space, as 
opposed to in a bubble wall. In addition, the final prediction for baryon 
asymmetry in the CKN scenario was very far from the equilibrium value (|^) 
because the sphaleron rate was slow on the time scales associated with the 
growth of the bubble. In our case, 6 changes slowly in time while the baryon 
number non- conservation is rapid. This allows a slow adiabatic adjustment 
of the baryon number to that which minimizes the free energy. 
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3 CP- Violating Operators and Coupling to 
the EW Gauge Fields 

In this section we describe some models of parametric-resonance-enhanced 
CP and B+L violation which are (at least immediately following inflation) 
spatially-homogeneous over each Hubble volume, but not fully-correlated 
over super-horizon regions. These models evade, in various ways, the usual 
problem that B-l-L generation in any one Hubble volume is completely un- 
correlated to neighboring Hubble volumes. In these models there is one or 
more feature which is totally-correlated (because of inflation), typically the 
infiaton VEV or the time-dependent part of Higgs VEVs as driven by cou- 
pling to the infiaton. But other features, such as initial values of other fields, 
may be totally uncorrelated from one Hubble volume to the next. Unlike the 
earher study [§] explicit CP violation is built in to the gauge-field equations 
of motion rather than just into initial values; this CP bias can overcome ran- 
domness due to initial values and to stochastic behavior of the solutions to 
the equations of motion. 

The CP-violating terms come from various higher- dimensional operators 
which can be generated from strong CP violation or from multi-Higgs mod- 
els showing spontaneous CP violation [0. The coupling strengths of such 
CP-violating operators need not be small, since in the early universe strong 
CP violation can be very much out of equilibrium. Generally the higher- 
dimension operators come from fermion loops; in the kinematic situation of 
the present paper, there is no way for a purely bosonic system of gauge fields, 
infiaton, and Higgs to lead to CP-violating terms in the gauge equations of 
motion. We will only consider local CP-violating terms, which would be ap- 
propriate at the beginning of pre-heating when the universe is very cold; in 
principle, modifications of locality can occur when the time dependence of 
various quantities begins to probe the fermion loops carrying the CP viola- 
tion. We will not consider that case here. 

There are two generic CP-violating operators which we will consider. The 
first is of the form 

kFTiG^^G"^ (10) 

where G is the dual field strength for the EW SU (2) groupQ, F{x) is a (possi- 
bly composite) gauge-singlet field and k a coupling strength with dimensions 

^We ignore hypercharge couplings in this paper. A factor of the gauge coupling g is 
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of mass to some negative power (depending on F). As usual, integration by 
parts of equation ([T0|) shows that the action can contribute to the gauge- 
field equations of motion only if F has a non-trivial dependence on x; we 
will seek for this dependence in out-of-equilibrium oscillations of the inflaton 
or other fields. We note that coupling of the form (|ToD as a source for topo- 
logical number density has been considered in a model with a cosmological 



pseudoscalar field coupled to hypercharge ||20 
Another possibility is: 

k' F Ti G ^yGyaGoi^i- (11) 

In this case, it is not required that there be any extra field F{x), since the 
operator in equation (|lT]) is not a total divergence. Nonetheless we include it, 
since there are terms in the gauge- field part of ( |lTl) which are total divergences 
and which could be important. 

When the action constructed from operators of the type ([10| , [Tl|) is added 



to the EW gauge action and the Higgs equation of motion is added, the 
dynamics studied in p are modified. We discuss the new equations of motion 
following some brief remarks on the physics behind the operators in (|TU|). As 
for the higher- dimension operator in equation (0), it is to be expected at 
some level whenever the lower-dimension operator in equation ( PUj) appears. 



4 CP- Violating Physics 

The CP-violating physics we are concerned with may come from out-of- 
equilibrium strong CP violation, from CP violation in the Higgs sector (with 
two or more Higgs fields), or from other causes. (Note that CKM phase 
effects are not strong enough to drive CP violation in the Standard Model.) 
We will discuss the first two explicitly. 

4.1 A model with strong CP violation 

Consider the CP- violating operator of equation (p!o|) , which has a typical 
axionic form, although we do not associate the field F with an axion. We 
assume that previous physics associated with earlier times has left a universe 

included in our definition of EW potential and field strength. 
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which has substantial strong CP violation in the QCD sector. As a proto- 
typical example, consider the Vj field, which is coupled to the QCD quark 
anomaly and to the gluonic topological charge density. This coupling to 
quarks is of the usual form 

M,gexp[^]g (12) 

(with suitable but irrelevant normalization of -Fr?')) showing that Vj jFr^i is an 
angle with period 27r. The potential energy of the r( field comes from (|1^) 
and from coupling to gluons; the latter must refiect the Witten-Veneziano 
[ plU relation and other requirements related to the 6'-angle dependence of 
QCD. The result is the standard form: 

V{tI) = mm, sqcd cos[^(^ + 27rj + V/F,,)] (13) 

where j = 0, 1, 2, A^c = 3 is the number of colors, 6 is the usual QCD vacuum 
angle, ^qcd is the QCD gluonic vacuum energy density (^qcd — 5's(TrG^) 
in terms of the strong coupling gs and field strength Gs] note that £:qcd < 
with our conventions using antihermitean gauge fields). To be specific, let 
us suppose that the vacuum angle is zero and that r]'/ F^/ is small enough so 
that we need consider only the j = sector explicitly. Now suppose that 
through the operation of some early-universe physics the field rj' deviates 
substantially from its equilibrium value of zero. Furthermore, this deviation, 
because of infiation, is (roughly) the same across the entire universe. Just as 
the infiaton does, the rj' field will begin to roll toward its equilibrium value 
(at which point strong CP violation is absent or small) on a QCD time scale 
~ GeV~^, a time scale long compared to all other time scales in the problem. 
The rj' field couples to the EW fields Gfj^i, through, e.g., quark loops, yielding 
an effective coupling 

^QCDr/'-^TrG^.G'^^ (14) 

It is straightforward to check that Kqcd order a QCD mass M of 1 GeV 
or so. The coupling to the EW potential requires an integration by parts, 
yielding a term in the action ~ kqcbV'^'^j where is a constituent of the 
EW gauge potential (see equation (^) below). 
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It might also happen that the field F coupled to the topological charge 
density oscillates at the infiaton rate, because of Higgs couplings to the infia- 
ton. For example, with the notation H for an electroweak Higgs field doublet, 
some unspecified physics may lead to an action of the type of equation (p!0|) 

2 

/ d^^H^HTiG,,,G^^ (15) 

where the field F = WH oscillates at infiaton rate scales because the Higgs 
field is coupled to the infiaton. 

4.2 CP violation in a multi-Higgs sector 

Since we are in no position to be very specific about CP violation in the early 
universe, we will illustrate the concept by using a model (T. D. Lee in |T0[|) 



which cannot account for all CP violation observed in the K — K system, but 
which can nevertheless represent an additional source of CP violation in the 
early universe. In this model there are two Higgs doublets Hi,H2, and the 
potential is chosen so that even though all couplings are real, the H2 VEV is 
complex and of the form 

m = ( ) (16) 

with real positive V2- The phase angle 6 arises from terms in the Higgs 
potential of the form 

V = H\H2 {DH\H2 + EH\Hi) + H.c. + ... (17) 

which can be rewritten to show terms involving cos 9: 

V = 2Dvlvl{co^e + ^^f + ... (18) 

Recall that the essence of EW-scale pre-heating is coupling of the infiaton 
to the EW Higgs fields. In this case, we invoke (for no deep physics reasons) 
a coupling of the form 

Ga^HlH2 + H.c. ~ ViV2a^ cos 6, (19) 

assuming for simplicity that the coupling G is real. Once pre-heating sets in 
and the infiaton a{t) begins oscillating, the angle 6 also begins to oscillate. 
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The CP violation is coupled through SU{2) fermions to the EW gauge 
fields, and one readily checks that the lowest-order fermion loop graph gives 
rise to a coupling of the type in equation (|TD|) which is proportional^ to 
V1V2 sin 6'. The oscillations of 6 then give rise to a non-trivial coupling to the 
EW gauge field 0. 



5 CP- Violating Equations of Motion 

There are several possibilities for CP-violating terms in the equations of 
motion, depending, for example, on whether the field F{x) in equation (p!0|) 
oscillates at the infiaton oscillation frequency, or at a much lower frequency. 
Before writing down these terms we set the notation by giving the Standard 
Model action for the gauge and Higgs fields, including a coupling of the Higgs 
field to the infiaton field a (but not writing other infiaton terms): 



S= (Tx 



^ Tr G^^G^'^ + D^H^D^H - X{H^H - K)^ f - G^a^H^H 

(20) 

We use, as before [|], the following ansatz for the EW gauge potential, which 
is g times the canonical potential and is written in antihermitean matrix form 
(so that D^ = + A^): 

Ao = 0; A, = (g)0(t), (21) 
with corresponding electric and magnetic fields: 

= Go, = (|)0; B, ^ \t,uGu = (g)0^ (22) 
For the Higgs field H we use the Standard Model form|] 

m = ( ) (23) 



V 



^The Higgs fields change left fermions to right, so both Hi and H2 must act. 

^The form of equation ( p3|) does not contain Goldstone modes, which contain vital 
topological information during actual creation of baryons; we will discuss these effects in 
a later section. 
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Evidently the topological charge density goes like 00^, which is, as it 
must be, a total time derivative. The topological charge density is the time 
derivative of the Chern-Simons density n^g. 

= &■ (24) 

The gauge and Higgs equations will be written in non-dimensional form, 
where the time variable is replaced by mt, the gauge field is replaced 
by 0/m, and the Higgs field v is replaced by v/vo- Here m = vog/2 is 
the vacuum W-boson mass in terms of the standard Higgs VEV Vq. The 
resultingi equations are: 



+ 20'^ + v^(j) + [CP - violating term] = (25) 



v + ^(l)\ + ip^)viv^ -l) + i^^)v = 0. (26) 



and 



where Mh = fov2A is the Higgs mass. An important property of equation 
( |26D is that if large (order one) values of are reached by resonant ampli- 
fication, the Higgs potential is strongly modified and v can approach the 
symmetry-restoring value of zero. 

There are, as discussed above, several possibilities for the CP-violating 
term. In these equations we will assume that the inflaton is oscillating sinu- 
soidally at a physical frequency u, corresponding to a dimensionless frequency 
r/2: 

rt, rm 
o- = o-Q cos(— ); uj = —. (27) 

Previously instead of using the Higgs equation of motion (^) the Higgs 
potential was set to zero, on the grounds that inflation stops, preheating 
starts, and the EW sector begins to undergo spontaneous symmetry breaking 
when the inflaton field reaches the critical value Oc- 



(28) 



^In these equations we ignore Hubble damping, which is miniscule, and damping by 
decay into fermions, which may have some importance at long time scales where other 
effects we ignore could also be important. 
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This led to 



V 



2 



0"n 



l + 2e + 2ecosrt; e= — ^. (29) 



Under the assumption that e is small ,| so that it can be ignored compared 
to unity in the constant term in (pQ]), the gauge equation of motion became 



+ 20^ + (1 + e cos rt)0 = 0. (30) 

Note that there is no CP violation in equation (^), since for any solution 
0(t) there is another solution — 0(t). CP violation can only come from initial 
values which favor one sign or the other of the CP-odd field 0. The modifica- 
tions we address in the next subsection add terms even in to the equations 
of motion, which therefore will have a built-in CP bias. And when we come 
to numerical analysis in Sec. 5, we will use equations (p5| , pB]) instead of the 
simplified form 1 + e cos rt in equation (|30D above. 



5.1 CP violation at the inflaton oscillation rate 

Before studying the complications of the combined Higgs-gauge equations 
|2^ ) which can really only be handled numerically, we will look at adding 
CP- violating terms to the simplified gauge dynamics of equation (|30D. This 
allows a certain amount of approximate analysis, similar to that used in the 
Mathieu equation, to be done. 

Using the Higgs dependence on time as shown in (^9|) and integrating 
by parts in the action of equation (|T5|) gives rise to an action of the form 
/ 0^ sin cut. We add the appropriate contribution to equation (|30D leading to 
the modified (and explicitly CP-violating) equation: 

+ 20^ + sin rt + (1 + e cos rt)0 = 0. (31) 

In this equation the parameter 5 is proportional to K^euv^. Note that if (p{t) 
is the solution to (0), then —(pit) is the solution to this equation when 5 is 
changed in sign. 

Without the 50^ term the equation of motion is essentially a Lame equa- 
tion, as discussed in [§. But with this term added we know of no way of 

^It is useful for approximate analysis to take the parameter e to be small, but it can 
be of order unity in actuality. 
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reducing the equation of motion to explicitly-soluble form. The major feature 
of the modified equation can be found by a standard Mathieu-like analysis, 
assuming that has the form 

(f){t) = a{t) cos(rt/2) + b{t) sin(rt/2) + c(t) (32) 

and assuming that a, b, c vary slowly on the time scale of oscillations. The 
new feature here is the presence of the secular term c, which is unnecessary 
when the term is missing in (0) . By ignoring all terms oscillating faster 
than the half-frequency r/2 one then finds: 

= rd + b[^ + ^-l-^ia^ + b^)]-5ac; (33) 

= -rb + a[^-'--l-^ia' + b')]~Sbc; 
—Sab 



l + 3(a2 + 62)- 

Note that since both 6 and c are proportional to the parametric-resonance 
parameter e, the new (last terms on the right-hand side of (^)) contributions 
are of order e^. The formal analysis is done in powers of e, which we therefore 
assume to be small. Then in (^) the equations for a and b are, to order e, 
the same as without the S term, and the secular term c is driven by the 
unperturbed values of a, b. These show standard Mathieu behavior of the 
form exp fit with 

/. = ±l[e^-(r^- 4)^/2 (34) 

(where we ignore temporarily non-linearities in the equations). Since e is 
small, the value of r must be nearly two for /i to be real and resonant am- 
plification to take place, which is to say that scales linearly with e. 

In studying the original equation with no 6 term but retaining the cubic 
non-linearity [§], it was shown that the quantity a'^+b'^ obeyed a certain equa- 
tion which revealed the conditions under which growth rather than damping 
(negative fi) occurred. With the 6 term added, this equation is: 

d ( 2 , .2-, -2ea& 25'^ab{a'^ - b^) 
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The second term on the right of ( pS] ) is of higher order and can be neglected, 
and then the conclusion is the same as before: Growth is only possibly if 
ah is positive (for negative e); otherwise there is damping. (When a, 6 are 
small, their growth leads to growth of the secular term c, at twice the rate 
of growth of a or 6; see (|33|) .) 
Writing 

a = y4cos\I^,6 = -y4sin\E' (36) 
with positive amplitude A, Ref. [|| showed from (|35|) that 

A = Aq exp ^* dt'-^ sin 2^(t') (37) 

with an accompanying equation for the angle ^ showing that it varied on 
the e rate scale, that is, slowly. Eventually ^ becomes large enough so that 
the product ah changes sign and growth turns into damping, or vice versa. 

What is happening to the Chern-Simons number in this model? The 
Chern-Simons number density is just 0^/87r^, and in the case where the 
CP- violating term proportional to 5 is absent, as is the secular term c, this 
quantity is purely oscillatory and has no appreciable long-term average or 
preferred sign. But things are different with 5 7^ 0; the secular average of the 
Chern-Simons density is: 

) = TT^cA^ (38) 



'M^' 167r 

Note that the sign of the Chern-Simons density is controlled by the sign 
of the secular term c. Under the usual assumption that the inflaton (and 
therefore the Higgs) field is correlated across the entire universe, because of 
inflation, the sign of 5 is also correlated across the universe and by (|33| ) the 
sign of c depends on the sign of the product ah, which may differ in each 
Hubble volume. But by equation ( P5| ) growth only occurs when this product 
is positive (since e is negative), even though a, 6 may separately be random 
in each Hubble volume. Of course, only those Hubble volumes where growth 
actually takes place have appreciable Chern-Simons density and in each of 
these volumes the sign of the Chern-Simons density is the same. There are 
Hubble volumes with Chern-Simons density of the opposite sign, but in them 
there is no growth and this density is small. The result is that at least for 
a while the entire universe has an average Chern-Simons density of fixed 
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sign, the same in each Hubble volume, even though the initial values of gauge 

potentials may be random in each Hubble volume. 

Numerical simulations discussed below show that on time scales of order 
the secular average can change sign, since on such time scales the product 

ab can also change sign. This may lead to relatively small Chern-Simons 

condensate values even in strong resonance. 

We cannot yet say how small or large this condensate value must be, in 

order to reproduce the observed value of B+L, or of the baryon-photon ratio. 

The smallness of the baryon-photon ratio must come, in our models, from 

a combination of the smallness of CP-violating parameters (such as occur 

in Higgs potentials), dynamical effects during pre- heating such as mentioned 

above, and sphaleron washout after pre-heating. 



5.2 Equations of motion with strong CP violation 

Since the 1]' field potential in equation (0) is determined by QCD parameters 
it is rolling very slowly on EW time scales and we can replace the time 
derivative of rj' by a term of order Mrj', where M is a typical QCD mass, and 
ultimately (after rescaling time and the EW potential (p as before) come to 
another CP-violating (p equation of the form: 

+ 20^ + 5'(p^ + (1 + e cos rt)0 = 0. (39) 

Here the parameter 6' is of order M/ m, the ratio of QCD and EW scales, or 
perhaps 10~^. 

As before, we must add a secular term c to the usual Mathieu ansatz of 
equation (|32D . Equations (^) are changed; the equations of motion for a,b 
are the same as with no CP violation (that is, set 5 = in the first two 
equations of (0)) while the equation for c becomes: 

- -^'C"^-^"^) (40) 



2[l + 3(o2 + 62)]' 



Provided that (as we have already assumed) 6' is roughly the same in all 
Hubble volumes, so is the sign of c, and so the Chern-Simons numbers in all 
Hubble volumes add with the same sign. 
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5.3 Higher- derivative CP violation 

We briefly consider here the consequences of adding the higher-derivative CP 
violation of equation (|lID . It is readily checked that this term yields actions 
~ / and ~ the latter is integrated by parts to ~ Hcjy'. So to 

the original equations of motion we must add terms of the form 

Fcj)^- -3(F02 + 2F00). (41) 

These discussed in connection with the numerical studies reported 

below, have dramatic effects because of the appearance of derivatives of (f). 

6 Numerical studies 

The range of parameter space is vast, and cannot be covered in depth here. 
We give several examples illustrating most of the effects discussed above. 
In these examples, we will use the following values of parameters unless 
otherwise specified: 

1. Gauge initial values: 0(0) = 0.01; 0(0) = 0.01 

2. Higgs initial values: f (0) = 0.5; ?)(0) = 0.5 

3. Epsilon parameter^ e = 0.5 

4. M^/2m^ = 1.5 (corresponding to a Higgs mass of about 140 GeV) 

5. Resonance parameter r = 2.2 

Other parameters will be specified as needed. Note that for this standard set 
of parameters the Higgs field is far from equilibrium; below we give numerics 
for initial values starting rather close to equilibrium. 

""^^We choose e positive for the simulations; this makes no difference to the ultimate 
interpretation of the numbers. 
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6.1 Sinusoidal CP violation with gauge backreaction 

Gauge backreaction is described by equations ( p5| , |26| ), to which we will 



add appropriate CP-violating terms. One of the most important features of 
this gauge backreaction is that it can help facilitate rapid Higgs transition 
through V = 0, which removes the sphaleron barrier and makes baryons from 
Chern-Simons number. The equations to be studied numerically here are 
(p6D plus a gauge equation with a sinusoidal CP-violating term: 

+ 20^ + v'^cp + sin rt = 0. (42) 

The initial conditions are as given at the beginning of this Section, except 
that 0(0) = —0.01; also, 6 = 0.51. This value of 6 is not necessarily realistic, 
but in general CP-violating quantities scale linearly in 6. The change of 
sign is done in order to get on the positive-e growth curve of the parametric- 
resonance instability of the simplified version of the gauge equation of motion, 
given in equation (PT|). Fig. |I| shows the results for 0, v, Chern-Simons density 
n^g, and the running average {n^^) of the Chern-Simons density, defined as: 

Ks) = 7 fdt'n^sit')- (43) 
t Jo 

There is no particular physical significance to the running average, but it is 
more convenient to display than the time-integrated Chern-Simons number, 
which covers too large a range to be displayed legibly. More to the point 
is a simple time integral of n^^ such as emerges from the dynamics, as in 
equation ( |60D below; this is 600 times the running average at the end of the 
time period simulated. 

One notes that in the present model the running average (n^^) itself 
changes sign from time to time on a long time scale, as would be expected 



from the analysis of Sec. |5.1| . So this model is not a particularly efficient 
way of generating B+L even in strong resonance. On the other hand, just 
how large a long-term Chern-Simons condensate needs to be to generate the 
observed B-l-L is not known, so this may or may not be a drawback. 

6.2 Constant CP violation with gauge backreaction 

In the present example we modify equation (p5|) by adding a quadratic CP- 
violating term with constant coefficient, as discussed in Sec. 



+ 20^ + v^cp + W02 = (44) 
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Figure 1: These figures siiow tiie beliavior of tlie dimensionless gauge potential 



Higgs VEV V, Chern-Simons density n 



, , and tlie running average of the 



Chern-Simons density (n) = (n^g) for the conditions of Sec. 3.1. 



We will take w = 0.5 for the numerical example. This is unrealistically 
large, but it makes it easier to see what is happening. Generally speaking, 
the long-term average Chern-Simons density is proportional to w; for exam- 
ple, as in the simplified analysis (with no gauge back-reaction) leading to 
equations which gives: 

or (n^s) 0.0008 for Figs. and g 

Fig. shows the time history of 0, v, the Chern-Simons density n^^, and 
the running average {n^g ) for initial Higgs values which are far from equilib- 
rium. The integrated Chern-Simons number is of order 0.1 for a CP-violating 
strength of order unity, and so might be of order 10~^ for more realistic CP- 
violating amplitudes w of order 10~^. This leaves room for inefficiencies in 
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Figure 2: These figures sliow tlie behavior of v, n^^g, and (n^g) for the 
conditions of Sec. with v{0) = v{0)=0.5. 



generating B+L from Chern-Simons number, washout of B+L, and other 
dynamical effects to reduce the baryon number to its actual level of perhaps 
10^^°. And, of course, there is no reason to believe that the numerical param- 
eters we use here apply to the real early universe. Experience with running 
many simulations of the type in Fig. |^ shows that the initial values of Higgs 
and gauge fields are not so important, and that final values of the running 
Chern-Simons average for unit-strength CP violation range from about 10~^ 
(somewhat smaller than the maximum expected, which is of order I/Stt^), 
to the non-resonant value of 0(O)^/87r^ which is of order 10~* for our initial 
conditions. Of course, these values would be reduced further by the strength 
of the CP-violating term, which is not small in our simulations. 

In the examples so far, the Higgs field has been far from equilibrium ini- 
tially and so it swings through the origin frequently. The next example, also 
for a constant CP-violating term, shows the importance of gauge backreac- 
tion in causing the Higgs field to oscillate through zero. For this example we 
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t t 

Figure 3: These figures show the behavior of (j), v, n^^, and {n^g) for the 
conditions of Sec. |6^ , with t;(0)=1.06, v{0)=0.1. 

choose Higgs initial values v{0) = 1.06, t'(O) = 0.1 with all other parame- 
ters as for Fig. The results are shown in Fig. |^. Because the Higgs field 
is near equilibrium, it stays for a while near unity, and the gauge field ap- 
pears to be trivially small. Actually (j) is slowly growing and finally becomes 
large around a time of 200. Fast growth of towards order unity values 
begins to send the Higgs field into wild oscillation, accompanied by growth 
of the Chern-Simons average value. We will see later that having the Higgs 
field go through the origin may be important for creating baryons from the 
Chern-Simons condensate. 
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Figure 4: These figures show v, n^^g, and {n^g) for the conditions of Sec. |6.3| . 

6.3 Higher-derivative CP violation 

Adding terms such as those in equation (^Tj) can and sometimes does lead 
to singularities in finite time because of the higher derivatives. We show one 
case where the singularity does not develop before dimensionless time 600. 
The equations to be solved numerically are the Higgs equation (^) and: 

+ 20^ + t;20 + w(a0^ + 00) = 0. (46) 

The coefficients w, a are -0.2 and 0.1 respectively. If w gets even a little larger 
in magnitude (say, w = —0.21) a singularity seems to develop, although we 
have not explored this phenomenon in any detail. Fig. ^ shows the results. 
The 0^ term in (^) is always of one sign and leads to secular CP violation; if 
only the (fxp term were kept, the running-average CP violation slowly wanders 
from one sign to the other. 
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7 From Chern-Simons Number to Baryons 



In this section we comment on some very difficult questions: How is the ho- 
mogeneous condensate of Chern-Simons number, |^ equivalent to a B-l-L con- 
densate, turned into actual baryons and leptons, which are localized states? 
What fraction of the condensate is turned into baryons and leptons? How 
does the winding number of the Higgs field enter in? We are in no position to 
give definitive answers; all these questions are still under active investigation. 

Let us begin with a comment concerning washout of baryons and lep- 
tons after reheating. The attractive feature of EW preheating is that the 
temperature after full reheating is rather smaller than the EW cross-over 
temperature T^, so that baryons and leptons created during preheating are 
not completely washed out (the sphaleron transition rate ~ exp(— Mg/T) is 
very small because the sphaleron mass Ms is large compared to the tempera- 
ture). We will assume that this mechanism protects baryons, once produced, 
from washout during reheating. Nonetheless, there can be washout during 
preheating, and we will have to discuss that. 

Our discussion will invoke insights derived from numerical and analytical 
work on these issues in the 1-1-1- dimensional Abelian Higgs model; from 
simple arguments based on effective temperature; and on approximations to 
an analysis based on modifications of earlier work on B-fL violation in two- 
particle collisions (where B-l-L violation is still extremely small because of 
the poor overlap of states, even though there may be no tunneling barrier). 

A real baryon is a spatially-localized state, so the transformation to 
baryons requires quasi-localized states resembling (but not necessarily iden- 
tical to) sphalerons; these states involve not only the gauge fields but also 
Higgs fields, which carry a winding number A^wind of their own. In the ansatz 
we have used so far the Higgs field has no winding number. 

Presumably a major influence on the dynamics of spatial localization is 
the instability of the gauge-Higgs equations to growth of spatial ripples, as 
alluded to in [^]. We will not discuss this mechanism here, which is non-linear 

Throughout this section we wiU make a distinction between total Chern-Simons num- 
ber 7V(-.g or Higgs winding number A'wind in a fixed volume V of space (or length L in 1-1-1 
dimensions) and their densities n^g, Jiwind, and similarly between the total energy EinV 
and its density e. Note that E/V is a function of n^g, riwind but not of V; nonetheless, it 
is appropriate to discuss the dependence of E on the total topological numbers N^^ , iVwind 
at fixed V. 
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and time-dependent (so it may be amplified by parametric resonance). 

7.1 Dependence of energy on topological quantum num- 
bers 

In the presence of a Cliern-Simons condensate it is energetically favorable to 
shift the winding number of the Higgs field Ny^md in the direction of Chern- 
Simons number. One way to see this is to note that the energy E = J d^xe 
depends on both A^'^.g and A^wind- In the homogeneous ansatz used so far 
the energy, with terms in 0^ and 0^, grows monotonically with n^^ since 
"■cs ~ The Higgs winding number iV^ind is zero (see equation (P^D). Now 
the energy of a state with finite values of A^wind, -^cs under a large gauge 
transformation with Chern-Simons number — A^wind, equivalent to one with 
zero Higgs winding: 

E{N^^ , AT^ind) = ^( A^cs - ^wind, ^wi,d = 0) (47) 
By the above remarks the energy with A^wind = has its minimum at A^^^g — 

N^ind = 0. 

More explicitly, once the system has become a condensate of more or 
less localized objects, these are described by a large gauge transformation at 
spatial infinity: 

A„ ^ Ud^U-\ (48) 

H(.)^±U[1). (49) 

The same gauge transformation U appears in both the gauge potential and 
in the Higgs field in order that the Higgs contribution / |-D0p to the action 
be finite in infinite volume.0 Evidently these localized states have lower 
energy than a homogeneous condensate which, so to speak, fills in the volume 
between the localized states. 

The vacuum states with A^^^g = A^wind are minima on the energy land- 
scape, so E{N^g, N^i^d) is zero along the line A^^^g = A^wind- Any transitions 

^^The gauge transformation U has the usual properties of large gauge transformations, 
including compactness on the sphere at infinity; a typical example is given in equation 
(H) below. 
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proceeding toward this line are energetically favorable. For example, begin- 
ning from the homogeneous ansatz state with large positive N^^ but zero 
^wind, it is favorable to decrease N^^ and increase A^wind- Of course, it is 
also energetically favorable to decrease N^^ toward zero with no change in 
-^wind) but then no baryons will be produced; B+L will rise and fall in lock- 
step with iVj^g, according to the usual anomaly relation, yielding no baryons 
at the end. So among the changes which decrease N^^ there is a competition 
between those which increase iV^ind and those which do not change A^wind- 

We illustrate these statements in Figs. ||, ||, ^ H. Fig. || is drawn 
assuming that the gauge boson mass (and a posteriori the sphaleron mass) 
are somewhere near their conventional values, so that there is a sphaleron 
barrier hindering topological charge change as indicated by the corrugations 
on the figure. However, as we have seen in the numerical studies, it is possible 
for this barrier to vanish as the Higgs VEV passes through zero. Fig. ^ is 
a view from above of the energy profile plot of Fig. ^. In Fig. ^ the filled 
circles along the N^^ = A^wind line indicate states of zero energy (with A^wind 
an integer), and the open circles on the dashed lines indicate no n- vacuum 
states of positive energy, constant along each dashed line. A conventional 
sphaleron transition from the point A moves along such a line with equal 
probability in either direction. This is not the case for topological transitions 
along other possible paths, such as the maximal energy gradient path (on 
Fig.|^, this path is orthogonal to the conventional sphaleron-transition path). 
Fig. 0, showing two cross-sections of the energy profile at A^wind = 0, 1, 
illustrates some possible transitions which change the topological numbers 
and the associated changes in energy. Fig. |] shows the energy profile at 
constant A^^s' with the corrugations corresponding to tunneling barriers. 

The energy profile E{N^g, A^„ind = 0) can be estimated in the case of 
static or slowly- changing spatially- homogeneous fields (see equation (^). For 



the 1+1 Abelian Higgs model |22| one has: 



EM/L = lvHgA,r = ln^^nl (50) 



and for the present 3+1 case 

EM/V = ^ {{8n^n,,r/' + (SttVs)^/^) • (51) 
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Figure 5: The static system energy E as a function of the Chern-Simons number 
and the winding number A^wind- Note that E{N^g, N^iad) is minimal along 
the vacuum line A^'^^g = A^wind and increases with | A^^^^g — A''„ind | • 

This profile is repeated along the A^wind axis with a corresponding shift along 
iV(^g axis so that the minimum i^^ = is always on the N^^ = iV^ind hne (see 
Figs. §,§). 

As mentioned before, there is no energy barrier for movement along the 
N^g axis, so this movement is completely unsuppressed and controlled by the 
dynamics of the gauge field condensate according to equation (0). This 
lack of suppression means that any (potentially large) number of baryons 
produced by variations in N^^ would instantly disappear immediately after 
returned to the energetically favorable value N^^ = A^wind, provided 
there were no transitions in Higgs sector (A^wmd remained unchanged). In 
other words, even though the fermions are coupled through the anomaly 
to Chern-Simons number, long-living fermionic states appear only due to 
transitions in Higgs sector which modify the energy profile along A^^^^g axis by 
moving the center (minimal energy) value of A^^^g around which it oscillates. 



29 




Figure 6: Possible patlis of topological transitions are shown on the A'^^.g - A^wind 
diagram, viewed from above. Path AB represents the simplest non-sphaleronic 
topological transition in the Higgs sector. See the text for explanation of the 
notations. 

As opposed to the absence of barriers for variations in A^^^g only, tran- 
sitions along the iV^ind axis require passing over sphaleron-like barriers (for 
simplicity, the barriers are shown at half- integer values in Fig. Precisely 
how high this barrier is depends on several circumstances, to be discussed 
below, notably the Higgs VEV and the spatial scales of the gauge and Higgs 
fields. In 3+1 dimensions spatially-homogeneous fields have a high barrier 
(provided the Higgs VEV is not zero), by an argument similar to that given 
earlier: it would be energetically favorable to separate regions with nontriv- 
ial gauge and Higgs fields in sphaleron configuration). In other words, the 
sphaleron, of size ~ m~^, corresponds to the minimum barrier height and 
the barrier height is greater for all other spatial scales. However, in the 
1+1 Abelian Higgs model the barrier at any N^^, N^ind is the same as the 
sphaleron barrier, as can be shown by making a large gauge transformation 
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Figure 7: A cross-section of the energy profile at constant A^wind- The lack 
of any barriers makes the downward evolution of N^^ completely unsuppressed. 
States A and B correspond to the same states in Fig. ^. At non-zero A^(,g the 
transitions which decrease {N^^ — A"„ind| become energetically favorable; a simple 
case discussed in the text corresponds to the path AB with N^^^ held constant. 



that eliminates the background gauge field and taking into account that the 
gauge field component of the 1+1-dimensional sphaleron is inversely pro- 
portional to the spatial volume. So the contribution of the gauge field to the 
(static) energy vanishes at infinite volume, while the Higgs field contribu- 
tion is identical to that of the sphaleron. The 3+1 dimensional case will be 
discussed below. 

In the homogeneous case discussed in previous sections, when the dimen- 
sionless gauge potential is of order one there is both a large Chern-Simons 
number and a large energy, but zero Higgs winding number A^'^ind- It is possi- 
ble that the energy in a volume containing one unit of Chern-Simons number 
can exceed the usual sphaleron energy, and (as we show in the next section) 
this does mean that transitions which change the Higgs winding number 
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Figure 8: Energy profile E{N^g, Ny^ij^^) along the AT^ind axis for several values 
of A^Qg. Note that iV„ind is integer for any field configuration; the curve profile at 
non-integer iV^ind reflects the presence of sphaleron-like energy barriers separating 
the states with AiV„ind = ±1. These barriers disappear at sufficiently high \N^g — 

A^wind I ■ 

become unsupprcssed. Note also that literal sphaleron transitions, which 
change iV^-,g and A^wind equally, on average do not change the Higgs wind- 
ing number, which is what is required to produce baryons. These do not 
change E{N^g — iV^md) and so their rate is not aflPected by the presence of a 
nonzero N^^ density. At any A^^g, the rates of transitions with AA^'^.g = -|-1 
and AA^pg — —1 equal each other and the contribution of normal topological 
transitions to total baryoproduction is zero. 

7.2 Effective-temperature estimates of transition rates 

Here we give some simple estimates of transition rates based on a commonly- 
used approximation, in which the net rate of baryon production is evaluated 



32 



from a first-order expansion of the rate in powers of an energy difference, 
leading to a form in which the rate is evaluated as a zeroth-order sphaleron- 
like part times a certain energy difference. In order that this be useful, one 
requires explicit knowledge of this energy dependence. A further simplifi- 
cation, used here, is to assume a Boltzmann dependence with an effective 



temperature. In Section ^]5| we give an approximate study of the sphaleron- 
like part of the rate, looking for conditions under which the tunneling barrier 
to baryoproduction can vanish. 

To evaluate the transition rate in the Higgs sector, assume that the tran- 
sition from zero winding to A^wind = 1 proceeds along Awind axis on the path 
AB (see Fig. Then the barrier height along this path can be estimated 
as 

= E,ph - (EiN^,) - EiN^, - 1/2)) (52) 

Similarly, the height of the barrier in the opposite direction is E_ = E^p^ + 
{E{N^g — 1/2) — E{N^g — 1)). If one could further assume a Boltzmann 
distribution with an effective temperature Tes, one would obtain different 
rates r± for transitions in the two directions, with a net rate: 

B/K^«^r,-r..Ifi-^ (53) 

where Fsph is the sphaleron rate at zero energy change. Then with the help 
of equations (|50| , |5lD , the winding number rate of change in 1+1 dimensions 
is: 

AT^a/L = Svr^r^p,^^ (54) 



and in (3+1) dimensions: 



iVwind/V = ^f^-^ (2(87r^n,3)V3 + M^(8vr^n,3)-V3) (55) 

Jeff Fes I '^'^W 

(here we assume the barrier height at A^ind = equals -E'sph, and note that 
the total rate is Fsph = F_,_ + F_). 

It is worth noting that the expression ( pSj ) becomes infinitely large at small 
n^g. This divergence appears because in 3+1 dimensions the homogeneous 
ansatz (|20|)-(|23|) used throughout the present paper is inadequate for field 
configurations with small topological charge densities, simply because gauge 
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field configurations with unit Chern-Simons number are localized in space 
(this is not the case in 1+1 dimensions). 

Strictly speaking, it is unclear whether the topological transitions in Higgs 
sector proceed along the A^wind axis or along the maximal gradient direction 
orthogonal to the iV^ind = -^cs However, the energy gains E± — E^ph 

in the latter case are larger only by factor of 2, which makes no qualitative 
difference in our analysis. 

With the increase in A^^^g — A^wind the energy gain obtained by transitions 
with A^wind = ±1 also increases, so at very large values of A'^^g — A^wind transi- 
tions along the A^wind axis also become unsuppressed because the energy gain 
( |52D becomes equal to or greater than the barrier height, i.e., the energy of 
a sphaleron-like Higgs configuration. The critical values can be estimated by 
putting E_ = in ([5^ ) and using equations (|50| , pID : 



~ = ^L.,,-' (56) 

In proper dimensionless units [E2] equation (|56|) becomes ti^'I** ~ 1/10 or 



Ai ~ 1/2) in 1+1 dimensions and 



< - (57) 

in 3+1 dimensions, which is about l/(27r^) per sphaleron volume. Equa- 
tion (0) is equivalent to 0>1, comparable to the maximum amplitudes of (j) 
seen in the numerical simulations. 



7.3 Baryoproduction 



Let us give some simple estimates of baryon production in the 1+1 Abelian 
Higgs model.Q Taking into account that the number of long-living fermionic 
states is equal to the net shift of Higgs winding number, and using equation 



41), one obtains 



B = 37T 



2 -^sph Tsph 
Mh Teff 



(58) 



The estimates in sections 7.3 and 7.4 can be straightforwardly generaUzed to the 3+1 
dimensional case using equation (psh, which should be applicable at least for large n^g. 
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and 

AB = 3n^^ f dt^n,,L (59) 

where AB is the total number of generated baryons. 

Ahhough the detailed discussion of Fsph and T^s time dependence goes 
beyond the scope of this paper, there are no reasons to expect any quali- 
tative difference from previous studies ^, |ll[. There it was argued that 



rsph/^eff is a smooth function of time with a sharp increase at the beginning 
of preheating and a slow decrease in the course of thermalization to the fi- 
nal reheating temperature. The exact timing depends on initial conditions. 
However, equation (|5^) provides a simple way to check how efficient the bary- 
oproduction is for any reasonable time evolution of Fsph and Teg. As long as 
the frequency of N^^ oscillations is much larger the inverse thermalization 
time, we may substitute Tsph/Tcs by 1 in the right-hand side of equation (^) 
which turns into the integrated average of Chern-Simons number density 

(n^,it))= f dtn^^it) (60) 







Time intervals when {n^^{t)) remains constant give zero baryoproduction; an 
increase or decrease in {n^^{t)) corresponds to a nonzero secular average of 
N^^ and thus to positive or negative baryoproduction (provided the transition 
rate rsph(t) is nonzero). The simplest way to get a nonzero A^'^^g secular 
average is to introduce a bias or tilt in E{N^^) through CP- violating terms 
either directly coupled to A^^^g in the bosonic Lagrangian, as in equation ([To|) , 
or coupled to the fermion density operator of equation (^ which modifies 
E{N^^) through the anomaly. It is also possible to shift the time average of 
N^g in other ways through dynamical effects already described in section ^. 



7.4 Washout from fermion backreaction 

Washout is one of the most prominent forms of fermionic backreaction. As 
long as the topological transitions keep going, newly-created baryons tend 
to disappear through the same anomaly mechanism that created them, to 
the extent that the decrease of fermionic density reduces the system energy. 
Washout (in 1+1 dimensions) can be accounted for by adding a dissipative 
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term to equation (|58D: 



B = 3n^^^n,,L-TnB (61) 

^ylH Jeff 



where the dissipation rate is generally proportional to the sphaleron rate 
Fsph. In 3+1 dimensions |Q, 



B 



39 Fsph 



(62) 



The solution of equation (|6TD can be found in the form (compare to equation 

my- 

AB = Sn^^e-'^'^ f dr e^^^^^-;^^n^,{T)L (63) 
Mh Jo reff(r) ^ ^ ^ 

(here 7(t) = /q dr Fb(t)). Direct use of equations ( |6T| , |63D appears to be prob- 
lematic in our simulations for two reasons. First, explicit time dependences 
of the rate Fgph and the effective temperature T^g are controlled by dynamics 
in the inflaton and Higgs sector which are beyond the scope of the present 
paper. Second, equation ( pT] ) assumes the fermions to be thermalized, which 
obviously is not the case if Chern-Simons number oscillates with frequency 
of order of m. 

However, it is easy to see that solution (^) decays exponentially with 
time if n^^ = and F^ 7^ 0. Therefore, it is desirable to keep baryopro- 
duction going until the topological transitions end (e.g. because of freeze- 
out during thermalization to a low reheating temperature). Otherwise, if 
baryoproduction ends before the sphaleron transition rate approaches zero, 
wash-out could crucially affect the final density of baryons. 

Again, the time dependence of the n^^ running average (|3) provides 
important information about the baryoproduction period and allows one to 
estimate acceptable thermalization times. For example, short-term baryopro- 
duction followed by stabilization of (^cs(^)) ^ certain value (see Fig. |^) 
would seldom survive the wash-out, while steady baryoproduction (such as 
is shown in Fig. ^ leaves ample time for the sphaleron transitions to freeze. 
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7.5 Dependence of tunneling barriers on and scale 
sizes 



Our purpose here is to find approximately the conditions under which the 
tunnehng barrier in the sphaleron rate T^ph = r+ + r_ (see equation (|53D ) 
can vanish. It will not be necessary to use an effective-temperature approx- 
imation. The techniques used here, based on work of a decade or more ago 
[P^ , PU| , P7| , [T^ on topological charge- changing transitions, are qualitative 
but useful. They go beyond 't Hooft's original work, whose famous tunnel- 
ing factor of exp[— 87r^/(7^] holds only at zero energy and, at that energy, is 
(classically) independent of scale size. We give an approximate barrier-factor 
formula with explicit dependences on energy, scale size, and Higgs VEV, and 
outline the regions of this parameter space where the tunneling exponent 
can vanish. We do not discuss the origin or distribution of spatial scales, 
but simply assumes that these arise through some mechanism such as the 
unstable growth of spatial ripples [If. 

There are other interesting approaches to this kind of problem, which 
we intend to investigate in the future. These include studies of non-Abelian 
gauge dynamics in the presence of a non-vanishing topological charge density 
or in the presence of a time-varying electric background potential 



and a study of the conversion of a time- varying background topological charge 
to fermions through the anomaly [Q. These works use various specialized 
backgrounds not exactly comparable to the Chern-Simons condensate used 
in the present paper, but they should still have qualitative applicabihty. 

There are at least two mechanisms which can remove tunneling barriers, 
and these may operate at the same time. The first mechanism is gener- 
ally important for large scale sizes (large means compared to the vacuum 
W-boson mass); it involves swinging of the Higgs VEV through zero. The 
second mechanism involves selection of a spatial scale such that the origi- 
nal Chern-Simons condensate is well-matched to baryon production through 
the anomaly. As one might expect, the best overlap occurs when the size 
scale is about m~^, and the corresponding configurations are somewhat like 
sphalerons in size, but it is apparently possible for tunneling barriers to be 
overcome at considerably larger spatial scales. However, baryon production 
is inefficient at these large scales. 

First we discuss the conditions for no tunneling barriers with a fixed and 
finite W-boson mass m, and then remark briefiy on what happens when there 
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is no mass because the Higgs VEV vanishes. 



7.5.1 The no-barrier condition for finite W-boson mass 

There must be a transformation of Chern-Simons number from the spatially- 
homogeneous condensate to one characterized by spatial inhomogeneities 
such as sphalerons. That is, in the expression for A^^^^g 

^cs = ^ / d^^ei.k Tr [AdjAk + '^AAjA^] (64) 

the relevant sphaleron-like configurations of the Ai are typified by gauge po- 
tentials which approach at infinity pure-gauge terms carrying winding num- 
ber; the Higgs field carries the same winding number for a minimum-energy 
configuration. The gauge in question was termed U in equations ^Qf ). 
A typical form for a gauge matrix U carrying unit winding number is 



(65) 



U = exp 



^T-rP{r,t)]; P{r,t) = 2 arctan[r/A(t) 



where X{t) is a monotonic function of t going to — oo at t = — oo and to oo at 
t = oo. A gauge matrix for Chern-Simons number A^^^^g is simply a product 
of N^g terms of the form (|65D translated to various spatial and temporal 
centers. We will refer to the spatially-homogeneous gauge potentials dealt 
with in earlier sections as form potentials, and the spatially-dependent 
sphaleron-like configurations as U form potentials. Of course, either form is 
at best an approximation; the (j) form potentials will develop spatial gradients 
by several mechanisms, and the U form potentials will not literally be of the 
approximate form we use below. 

Long ago, Bitar and Chang [^] constructed Minkowski-space gauge po- 



tentials of unbroken SU{2) gauge theory whose asymptotic behavior was 
precisely that of equation (|65D. They chose the single dynamical degree of 
freedom to be the function X(t) of this equation; the non-asymptotic gauge 
potential is further parameterized by a non-dynamical scale coordinate and 
translation coordinates. The parameterization is conveniently written as: 
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where p is the scale coordinate (we will not need to display the translation 
coordinates, which we set to zero). Bitar and Chang show that the 
topological charge of this potential is unity for any scale coordinate value. 
They also show that the dynamical degree of freedom A has a Hamiltonian 
which is quadratic in A but a complicated function of A, and computed the 
barrier exponent / pdq at zero energy from this Hamiltonian, recovering the 



't Hooft result. Later this work was extended |^ to EW theory with a 
Higgs field (as well as a chemical potential for Chern-Simons number); this 
work then formed the basis |]l3l for an investigation of scattering processes 
involving topological charge change at very high energies. In this paper, we 
extend these earlier works to cover the entire range of scale coordinates (Ref. 
[0 only covered the regime of small scale coordinates). 

We wish to find the conditions under which the corrugations or barrier 
factors described in connection with Fig. |^ vanish, within a framework gen- 
eral enough to go beyond thermal quasi-equilibrium. Begin with a general 
formula of Cline and Raby for the diffusive rate. Originally the Cline- 
Raby formula was given for thermal equilibrium conditions, but it is easily 
modified for non-equilibrium conditions. The derivation is slightly different 
from Cline and Raby's because of the non-equilibrium nature of the process. 
Consider diffusive dynamics of the Chern-Simons number N^^{t) in which 
the diffusion constant^ F can be written in the usual form: 

°° if 

(67) 

where V is the volume of space and the brackets refer to a trace over the 
density matrix; this density matrix is, for the present purposes, taken to be 
diagonal in the energies Ei,Ef of the states summed over, with entries p{i). 
In fact, the density matrix is changing in time and so does not commute with 
the Hamiltonian, but at the present level of (in) accuracy this is an inessential 
complication. Cline and Raby take the i states to be in states at t = —oo 
and the / states to be out states at t = +oo, and then use the formula: 

/oo 
dt{t,E,\N^^{t)\f,Ef) = 2n6{E,- Ej){t,E,\N^M\f,Ef) (68) 
-oo 



^''Roughly speaking the diffusion constant F used here is equivalent to ^Fsph of section 
71. 
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where Tji is the zero-temperature S-matrix element from the initial to the 
final state, and ANfi is the change in Chern-Simons number (or topological 
charge) from the initial to the final state. Substitution in ( pB| ) yields: 



'^4(p/-p.)(AiV^,)'l^./ 



(69) 



if 



Further progress depends on analysis of B+L- violating S-matrix elements 
at zero temperature, a subject of some considerable interest a decade ago 
(see, e.g., [0). Ref. |0 gives the following very crude approximation to 
the |AiV/j| = 1 S-matrix elements at fixed energy E: 



where Q is a barrier exponent, is the total number of particles involved 
in the scattering process, p is a scale collective coordinate, and kj is the 
magnitude of the three-momentum of particle j. This formula is based on 
a transcription of familiar Euclidean formulas for scattering in the presence 
of instantons to Minkowski space. Unlike naive instanton-based amplitudes, 
the above amplitudes Tij behave properly at high energy (where one expects 
p ~ N/E) but are not correctly unitarized; this will not be an important 
issue here.0 

The barrier exponent Q was originally |jl3| given for small p. The appro- 
priate expression for all p can straightforwardly found using the techniques 



F{o = ie + ir'/'ne + ir'^' + ^(t) w(e) - tt^^^p/i- (72) 



Here we make it explicit that the W-boson mass m{t) depends on time, 
because the Higgs field depends on time. The function /(^) is not expressible 
analytically, but one can show that / is positive and obeys /(^) < /(O) = 4/3. 

^^See [Sll for a multi-channel study of unitarization effects. 



(70) 



of H, HI: 




(71) 



where the function F is: 
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The barrier factor Q vanishes if the function F is always negative, which 
happens for certain regimes of energy E, mass and scale coordinate p. 

Clearly, there is always a finite barrier at p = (just the 't Hooft barrier if 
m=0). For non-zero p it is generally true that if = 0;p,m,E) vanishes 
then F < for all ^. The no-barrier condition then is: 

Am?p^ Epq'^ ^ 
1 + ^ < 0. (73) 

Consider first the usual case where m is the standard W mass. The minimum 
energy E yielding equality in equation ([73|) occurs at p = V3/2m, Emin = 
4n^y/3m/g^. This minimum should be the sphaleron energy, and it is indeed 
a very good numerical approximation ||26| in the limit of large Higgs mass. 



Next we take up the case where the W mass, or Higgs VEV, is near zero. 
7.5.2 Higgs VEV near zero 

It appears that large scales p, which would be expected in the first stages 
of transition of Chern-Simons number from form to U form, will have a 
disastrously large barrier factor going like exp(— const mp/g'^), if the mass 
m is anywhere near its standard value. However, if — as discussed in Sec. 
5A — the Higgs field v oscillates through zero because of gauge back reaction, 
the gauge mass m = gv/2 vanishes and the barrier Q will vanish periodically 
(see ( |72D at energies E ~ l/{pg^)] that is, more easily at large scales than at 
small.|3 So Higgs oscillations are one potentially-vital means of generating 
baryons from Chern-Simons number. It is easy to check that if the barrier 
factor is unity for a time r during a Higgs oscillation period Tniggs, the aver- 
aged barrier factor over a Higgs period is not exponentially small, but is of 
order r/rmggs. 

Unfortunately, this is not the end of the story, since there are other pos- 
sible suppression effects to deal with; a poor overlap of anomaly-produced 
baryons with the Chern-Simons condensate can be disastrous. 



^^The idea that preheating causes oscihations or vanishing of v and therefore reduction 
or ehmination of the barrier factor is given in Refs. M p2l. 
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7.6 Good overlap condition 

Here we explore how the conditions of both no barrier as well as a good over- 
lap between a U form potential and baryons produced through the anomaly 
can be satisfied. 

Even if v does not vanish during preheating it is possible in principle to 
have a zero barrier factor, depending on what values of the dimensionless 
gauge potential are reached. Values of |0| near or slightly greater than one 
are fairly readily gotten, as is evident from the structure of the equations of 
motion and from the numerical studies reported above, and for such values 
of (f) there need be no barrier. However, it is possible for |0| to be small 
compared to unity; we explore that possibility here. As one might expect, 
for values of rather less than one, it requires a rather large region to gather 
together enough energy to overcome the barrier. It will turn out that spatial 
scales much larger than the inverse of the vacuum W-boson mass are self- 
consistent only if the Higgs VEV does go near zero. 

Consider a region of space of size p as defined by gradients appearing 
during preheating and further unstable amplification. In this region the 
total energy and Chern-Simons number are approximately: 



Eip) ^ ^impm' + n NUp) - (74) 



Inserting these estimates in the no-barrier condition ( [73| ) yields: 

with the last form holding for small |0|. In this case the total energy and 
Chern-Simons number in terms of d) are: 



47r^ _7rV2 



E = —^: iV,„, = ^ = 0(1). (76) 



The appearance of inverse powers of suggests the inefficiencies of avoid- 
ing a barrier when parametric resonance amplification of the EW gauge po- 
tential is small: The energy is large compared to the sphaleron mass, but the 
change in B-l-L is only of order one (per unit fiavor). Still, these inefficien- 
cies could be tolerable in view of the exponentially-small efficiency of actual 
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barrier penetration. When amplification leads to |0| ~ 1, avoiding a barrier 
is rather like having energy at or above the sphaleron mass. 

Next turn to the conditions specifying a good overlap between the 
condensate and baryogenesis. We will find that a good overlap means |0| 
cannot be too far below unity, but we cannot quantify this statement with 
the present crude approximations. The point of a good overlap is that having 
an energy larger than the barrier energy is by no means sufficient in many 
cases to lead to unsuppressed B+L violation, as many authors have discussed 
1^, |13|. For example, in the formula ( [70|) for S-matrix elements, the other 



factors integrated over p in that equation can lead to suppression of the 
S-matrix element by a factor of order exp[— CSvr^/fyf^] with C, ~ 1/2 or so 
r^j. This sort of suppression even in the absence of a tunneling barrier 



comes about because of a very poor overlap between multi-particle initial 
and final states when these have very different numbers of particles and at 
least one particle number is very large. When the states are sufficiently 
similar there is no such extra suppression, which is what happens in thermal 
equilibrium at large enough temperature. In the present case something 



similar happens. The initial state in the Cline-Raby formula (|68| , |69D is not a 
conventional particle state; it is a coherent state somewhat similar to the 
form potential, but with a spatial size p coming from various effects, such as 
growth of unstable momentum modes . If this state has a large overlap with 
the U form potentials such as in (|66D, and if there is no barrier for tunneling, 
then the amplitude for baryon creation will be unsuppressed. One can say 
that the exponentially-small rate of B+L violation in 2 ^ collisions stems 



from the Drukier-Nussinov effect [0 that it is extremely unlikely for a two- 
particle collision state to couple well to a soliton like the sphaleron, but in 
our case the Chern-Simons condensate may, for certain values of and p,look 
enough like a "soliton" for there to be good overlap. 

We seek this substantial overlap between a generic form potential, some- 
how modified to have an overall spatial scale p, and a U form (Bitar-Chang) 
potential, when p is a large scale compared to other spatial scales. At large 
p the Bitar-Chang fields scale as: 

gE r^^- gB^ \. (77) 
p^ p^ 



43 



By comparisonQ^o the form potentials of equation (|2^ ) one sees from the 
B field consistency with the relation mp ~ 1/|0|, and that consistency with 
the E field can be achieved if A ~ In 0. 

Now return to the approximate form of the S-matrix elements (ffOl). The 
process to be described is the transformation of a form field with certain 
spatial scales to a final state of particles, including a set of particles which 
violates B+L. The process is only interesting if there is no tunneling bar- 
rier, so we assume the validity of equation (^). The next question to ask is 
whether one can, consistent with (|73D , argue for a non-suppressed overlap be- 
tween the modified form and the final state. In the crude approximation of 
equation (|70| ) this simply amounts to asking whether the constraint on scale 
coming from the final-state wave function is consistent with other informa- 
tion on that scale, such as the condition ([73|) for no barrier. In ([70D there 
occurs a product of exponentials of the type exp(— fcjp) where ki is the spatial 
momentum of the i^^ particle. It is reasonable to assume that all the particles 
are effectively massless, and then the product of wave functions reduces to 
exp(— ii^p). One easily sees that, if the barrier exponential factor exp —Q is 



unity, the p integrand maximizes at p = N/E, just as estimated |jT3[ for more 
conventional S-matrix elements. So a good overlap simply means that the 
number of produced particles is (more or less) determined. This turns out, 
taking into consideration the no-barrier condition (|76D, to yield ~ 1/(7^0^. 
(One might argue that in fact particles do have mass, and so one should 
have < E/M ~ l/g'^\(j)\ which would require |0| at least of order unity. 
However, m is the W-boson mass, and all baryons and leptons except for 
baryons with top quarks are not nearly that heavy.) 

As in [0 we can form a rate from the Cline-Raby formula ^ by multi- 
plying the squared S-matrix elements from ([70|) by massless A^-particle phase 
space [0. Nothing quantitative should be trusted about the resulting equa- 
tion (|78D below, except for its dependence on 0: 



V N\g^^ 'g 



2N ^n'iA.' ^(j)2' 



Here 0c is a critical value of separating small from large rates (in the 
present highly inaccurate approximation, 0^ = e~^). Whatever 0c is, and it 



^^The fields (j), (j) in equation ( P4 ) are dimensionful; the ones we use now are scaled by 
appropriate powers of m to be dimensionless. 
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is evidently of order unity, it is clear that for N >2 and for 0^ > 0^ the B+L 
process is unsuppressed, while in the opposite limit it is strongly suppressed. 
Note that if the dimensionless potential is of order unity the energy scale 
is the sphaleron mass and the spatial scale is m~^. So we conclude, not 
unexpectedly, that if the Higgs VEV is not oscillating near zero spatial scales 
near the sphaleron mass must be formed; if the Higgs field is oscillating near 
zero, much larger spatial scales will serve for unsuppressed B+L production. 

8 Conclusions 

In this work we studied a scenario of baryogenesis in the Standard Model, 
based on inflation on EW scales with the Higgs coupled to the oscillating 
inflaton in a preheating phase. If baryons could be generated, they could be 
saved from washout at reheating because the reheat temperature is less than 
the EW cross-over temperature where sphaleron washout can be large, so 
the topological transitions can completely freeze out before the preheating 
ends. However, the dynamics of conversion of Chern-Simons condensate to 
B+L and the dynamics of washout of this B+L are very complex and not 
yet well-understood. 

Although we have no rigorous proof that B+L production involving only 
standard-model fields (coupled to the inflaton) is truly viable, we can at least 
say that so far that we have not identified any mechanisms which rule out a 
pure EW/infiaton scenario. 

Extending the earlier work of Ref. P], we added explicit CP violation to 
the spatially-homogeneous classical gauge equations of motion, rather than 
(as in 1^) depending only on initial value of the functions 0,0 for CP vi- 
olation. The effects of explicit CP-violating terms was studied analytically 
in some simple models, as well as numerically. We studied numerically the 
classical dynamics of the Higgs field, including classical back reaction of the 
gauge field on the Higgs field. 

The gauge reaction on the Higgs equation of motion is important in two 
respects: It can broaden parametric resonances by getting away from pure 
sinusoidal variation of the Higgs field, which was the only case considered 
in 0. And it can, as discussed above, lead to conversion of Chern-Simons 
number to baryons which for some fraction of the time is unsuppressed by 
tunneling barriers, since gauge backreaction often modifies the Higgs poten- 
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tial and causes the Higgs field to oscillate through zero, rather than staying 
at the bottom of the potential well. 

There are various forms of CP violation which could be added to the 
equations of motion; we explored three in the gauge equations and one in the 
Higgs equation. In the gauge sector one came from strong out-of-equilibrium 
CP violation; another came from a multi-Higgs sector with spontaneous CP 
violation; and the third was a higher-derivative CP-violating term, which 
could be associated with other lower-dimension operators. In the Higgs sector 
we explored spontaneous baryogenesis in a multi-Higgs model, leading to an 
effective chemical potential for baryons. Generally speaking, the resultant 
secular average of the Chern-Simons number is a few orders of magnitude 
less than the dimensionless coefficient multiplying the CP- violating term in 
the gauge equations of motion. 

Formation of a spatially-homogeneous Chern-Simons condensate is only 
the beginning of the story; it is necessary to convert this condensate into 
something resembling a condensate of sphalerons of the same Chern-Simons 
number, in order to generate the Higgs winding number which is converted to 
baryons. The usual approach of invoking a thermal (or effectively thermal) 
regime immediately following the parametric resonance regime is not likely 
to be applicable. A quantitatively-adequate study of these non-thermal non- 
equilibrium processes remains to be done, but we have given useful criteria 
for evading the two possible process which lead to exponential suppression of 
baryogenesis. The first, of course, is the topological-charge tunneling barrier. 
The second is a poor overlap between initial and final states in a process even 
with no barrier. Based on earlier work which studied the (im)possibility of 
baryon production in accelerator collisions, we gave crude estimates of the 
parameter ranges which avoid both barriers and bad overlap, and pointed 
out the infiuence of the Higgs VEV passing near zero for these parameter 
ranges. The conclusion is that if the Higgs VEV stays near its vacuum value, 
spatial scales near the vacuum inverse W-boson mass and energies near the 
vacuum sphaleron mass can avoid both forms of suppression, quite unlike 
the 2 ^ N collision case where poor overlap cannot be avoided. If this VEV 
goes near zero, the spatial scales can be larger and the energy scales smaller. 

If this view proves to be correct, the rate-limiting step in baryogenesis 
will be the conversion of form EW fields to U form fields, by growth of 
perturbations at various spatial scales. 

The scenario of regions of almost spatially-homogeneous Chern-Simons 



46 



condensate of order of the EW Hubble size, may have consequences beyond 
B+L generation. It has been proposed |33| that any mechanism of B+L gen- 
eration involving the EW anomaly will leave its trace on the early universe 
through large-scale helicity of the Maxwell magnetic fields which descend 
from EW gauge fields as the universe cools. As is well-known, if these mag- 
netic fields have typical EW scales of 10~^^ cm, the present-day magnetic-field 
scale will be far too small even taking expansion of the universe into account. 
Following earlier studies in magnetohydrodynamics, Ref. |Q proposed that 
the helicity would drive an inverse cascade to longer scales. This idea was 
further pursued by Son and by Field and Carroll |Q. If the EW fields 
which lead to Maxwell magnetic fields are generated on the EW Hubble size 
scale, one gains many orders of magnitude toward seeding present-day inter- 
galactic magnetic fields by EW processes. Moreover, there might be much 
more Maxwell helicity than the minimum required to produce the present- 
day abundance of baryons, since only a fraction of the form Chern-Simons 
condensate will be turned into baryons. The Maxwell helicity produced dur- 
ing preheating amplification of CP violation can be very much greater than 
the number of baryons. Work is underway to investigate these points. 
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